Thermodynamic Compressibility of a Two-Dimensional Electron System: Signature of 

a Droplet State 
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We have used a field-penetration method to measure thermodynamic compressibility of a moder- 
ately interacting two-dimensional electron system (r^ « 0.5-3) in a three terminal GaAs/AlGaAs 
device, fabricated with an epitaxial lift-ofi^ technique. We found that the density and temperature 
dependencies of the compressibility are qualitatively different from that observed in earlier studies 
of the 2D hole system, where interaction energies are considerably stronger. We show that the 
observed characteristics can be described by the recently developed formalism for compressibility of 
the droplet state. 



The apparent metal-insulator transition discovered in 
low-density, high mobility Si, SiGe and GaAs based two- 
dimensional semiconductor systems continues to capture 
much interest in the general condensed matter physics 
community [1]. Although the experimental work in this 
research area has been predominately transport stud- 
ies, there is recent work, both experimental [2-4] and 
theoretical [5,6] devoted to the study of the thermody- 
namic compressibility of the interaction of disordered 
two-dimensional systems, in the context of the 2D metal 
- insulator transition. 

In our earlier measurement of a strongly corre- 
lated two-dimensional hole system (2DHS) in a p-type 
GaAs/AlGaAs heterostructure, we found that the neg- 
ative inverse-compressibility at low-densities reaches a 
minimum value at the apparent metal-insulator transi- 
tion and than increases dramatically with further de- 
creasing density. The results show a strong correlation 
between the qualitative change in thermodynamic com- 
pressibility and the transport signature of the 2D metal- 
insulator transition. In that experiment, the 2DHS has a 
rather strong correlation due to the heavy effective mass. 
The dimensionless parameter which measures the ratio 
of the correlation energy to the Fermi energy ranges be- 
tween 5 to 25 and the qualitative change in the compress- 
ibility was observed at rg—lZ. To understand whether 
this behavior is specific to strongly correlated systems 
and how it is linked to the 2D metal-insulator transition, 
we have studied another two-dimensional electron sys- 
tem (2DES) where interactions are much weaker, — 
0.5-3. There is no transition in the transport. We found 
that the compressibility behavior in the present system 
is qualitatively different from that observed in the 2DHS 
system. Its density dependence can be well described by 
the theory [6] of percolating droplet states. 

To measure the thermodynamic compressibility, we 
have used a field penetration technique, developed ini- 
tially by Eisenstein. [8] Unique to our setup, the 2DES is 
sandwiched between two metal electrodes (the front and 
back gate) as shown in the inset of Fig. 1. By using 
an epitaxial lift-off technique for removing the MBE sub- 
strate, we were able to place a back-gate at a distance 
of 1000 nm from the 2DES, while the front gate is 167 
nm away. The samples used in this experiment were cut 



from an n-type GaAs/Al^jGai-^^As wafer fabricated by 
molecular beam epitaxy. A 70 A Al.30Ga.70As undoped 
spacer was used with an atomic planar doping layer of Si 



just above with a concentration of 1.5 x 10 



The 



rest of the sample, designed specifically to do the sub- 
strate removal, was identical to the one used for the 2DHS 
compressibility study [2] . The mobility of the sample was 
186, 000 cm?' /Vs at a density of about n — Qx 10^^ cto~^ 
when the sample was unbiased. The wafer was carefully 
designed so that the mobility of the carriers was compa- 
rable to that in the 2DHS studied earlier. 
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FIG. 1. Penetration current components Ix and ly vs gate 
voltage for five temperatures at an excitation frequency of 
100 Hz: the temperatures are labeled from right to left as 
0.28 K, 0.56 K, 0.83 K, 1.0 K and 1.2 K respectively. Inset: 
the schematic of the experimental setup. 



To perform the experiment, we applied a 10 mV AC ex- 
citation voltage Vac to the front gate. A DC voltage Vg 
was superimposed to control the carrier density. With 
the 2DES grounded, the penetration current from the 
front gate to the back gate was detected by a lock-in 
amplifier at an excitation frequency of 100 Hz. The fre- 
quency was chosen in such a way as to maximize our 
signal-to-noise ratio while remaining in a frequency inde- 
pendent regime. 

As described before [2], by modeling the system as a 
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distributed circuit, both the quantum capacitance Cq and 
the resistance Rs of the channel for the 2D carriers could 
be extracted individually based on the measured values 
for the in-phase and 90° phase ly current components: 
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Rs. 



(1) 



(2) 



where m is the frequency of the excitation voltage, and 
Ci and C2 are the geometric capacitances between the 
front and back gates with respect to the 2D layer. 

The quantum capacitance per unit area Cq is related to 
the compressibility, k = (^)/"-^, by Cg = e2(^) , where 

is the chemical potential and n is the carrier density. 

Figure 1 shows both the in-phase Ix and 90° phase ly 
components of the penetrating current as a function of 
the gate voltage for five different temperatures ranging 
from 0.28 K to 1.2 K. The d/z/dn value, which is in- 
versely proportional to the compressibility, is extracted 
and plotted in Figure 2 as a function of the value in the 
bottom axis. At 100 Hz, d/i/dn is almost directly pro- 
portional to Ix ■ There are important differences between 
this curve and the same curve for the 2DHS. First, it 
is immediately apparent from the data in Figure 1 that 
the resistive component of the signal does not have a 
temperature independent crossing point like that for the 
2DHS. This should be expected from a sample with only 
weak to strong localization crossover behavior as shown 
in the transport measurements. The transport proper- 
ties of the 2DES can be understood in terms of the more 
conventional scaling theory. [7] Secondly, d^/dn for the 
2DES gas is positive everywhere. At high densities, this 
dependence is expected. As one can see from the quanti- 
tative agreement with the Hartree-Folk (HF) approxima- 
tion calculation for low values. [9]. Since the effective 
mass of electrons in GaAs is only 0.067 times the rest 
mass of an electron (roughly 5 times less than the effec- 
tive mass of the holes in GaAs), the maximum rg value 
of this particular system is only about 3 for the low- 
est density studied here. Positive compressibility means 
that the kinetic energy is always greater than the ex- 
change energy for the present system. Third, although 
there is still a minimum in d///dn, just like that observed 
in the 2DHS, the turn- around is not nearly as sharp as 
is the case for that system and occurs at a density of 
roughly [ts ~ 1.4). Finally, the temperature dependence 
of dfi/dn is rather strong, in contrast to what is observed 
in the 2DHS where dfj,/dn is temperature independent for 
a large range of temperatures from 0.3 K to 4.2 K. 

Over the years, it has been suggested that the inter- 
acting 2D electrons in the presence of disorder can be un- 
derstood by the droplet, or puddle, model. [11,12] In this 
model, the electrons are separated into "liquid" droplets 
with local density higher than the average density and 




FIG. 2. dfi/dn is plotted against rs for the data in Figure 
1. The solid line is calculated using HF approximation. 



"gas" islands with lower density. For a given disorder, 
due to the potential fluctuation of the donors, as the 
Fermi energy is decreased, the "gas" region increases. In 
this model, a cross-over from metallic behavior to insu- 
lating behavior should be seen around the percolation 
threshold. [13] 

Specifically, the thermodynamic compressibility has 
been recently calculated by Shi and Xic for an inter- 
acting system within the local density approximation 
(LDA). [6] The LDA result produced a "turn-around" 
in dfi/dn which appears to be similar to that observed in 
the 2DHS system. [2] Furthermore, the LDA calculation 
was in good agreement with numerical simulations using 
"liquid" droplets and "gas" islands. 

Following the work of Shi and Xie, we have calculated 
the compressibility numerically using LDA for the exper- 
imental parameters and found that the basic character- 
istics of the electron data can be well described by the 
formalism of the LDA. 

In the LDA calculation, the local density (?ie//) rather 
than the average density, is used to determine the total 
energy of the system. The local density is assumed to be 
rieff — n/f{n), where f{n) is the fractional area of the 
high density region. The chemical potential is therefore: 



H{n) = d{eo{neff)n)/dn, 



(3) 



where Eq is the energy density for the uniform electrons 
and is given by [9]: 
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where e is the dielectric constant of the host material 
and m* is the effective mass of the 2DES. 

The chemical potential and d^/dn are calculated nu- 
merically. As shown in Figure 3, the LDA approximation 
describes the data rather well in the high temperature 
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limit. To fit the data we use the following relation for 
/("): 



1 



(1 + (no/n)«) 



(5) 



no=1.6 X 10^^ cm^^ and a=2.0 have been used to pro- 
duce the theoretical curve in the figure. It turns out 
that the numerical fit is very sensitive to the value of 
a and uq. Small deviations of about 10 percent from 
no— 1.6 X 10^^ cm^^ or a— 2.0 dramatically effect the 
curve in such a way that we can no longer qualitatively fit 
the data. Incidentally, the value of a of 2.0 is very close 
to that of 2.3 used to fit the numerical simulation by Shi 
and Xie [6] for off-plane, un-correlated charge impurities 
and this parameter should depend on the form of the dis- 
order potential. We do not understand why the data at 
lower temperatures or lower densities deviate from the 
theoretical curve. It is possible that the temperature 
dependence of the size of the lower density islands is al- 
lowing the field to penetrate to the back gate more easily 
at the lowest temperatures. This could create a temper- 
ature dependence in f{n) in equation 5, which could ex- 
plain a temperature dependent dii/dn curve. However, it 
is unlikely that the temperature dependence would come 
about from a single temperature dependent parameter in 
the above equation. [16] 



mathematical singularity, which of course is inconsistent 
with that for percolation. As an example, we plot a curve 
of the LDA using f{n) = (1 — (no/p)^). A reasonable fit 



is obtained with no— 3.5 x 10 



and f3 = 0.06. This 



observation implies that the physics of the 2DHS which 
appears in the large regime, CANNOT be described 
by the droplet model. 

The necessity of a singularity in f further suggests that 
there is a phase transition in the 2DHS system. This 
observation further supports the notion that the non- 
analytical behavior of the thermodynamic compressibil- 
ity is consistent with the unusual signature of the metal- 
insulator transition observed in the transport for the 
same sample. Indeed there are theories describing the 
2D metal-insulator transition as a phase transition from 
a liquid to a Wigner glass state [14,15] or from a para- 
magnetic state to a spontaneously ordered ferromagnetic 
sate [17]. It is difficult to determine the existence of these 
potential states from these measurements due to the ab- 
sence of a quantitative analysis of the compressibility of 
the Wigner glass state or the ordered ferromagnetic state. 
For example, the compressibility of a Wigner crystal, in 
the absence of disorder, is known theoretically [18] to be 
only twice that for a liquid in the HF theory. However, 
compressibility for a glass state, in the presence of real- 
istic disorder, has not yet been calculated. 
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FIG. 3. The data is compared with the LDA approxima- 
tion. The LDA approximation describes the data well in the 



high temperature limit, 
used for the fitting. 
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FIG. 4. In contrast, the LDA approximation cannot de- 
scribe the dfi/dp vs. p curve obtained in the 2D hole system, 
since the theoretical curve varies slowly with p. A reasonable 
fit can only be obtained if there is a singularity in the filling 
factor /(n), as described in the text. 



It is important to note that the accurate data fit in 
Figure 3 to the theoretical curve should not give the im- 
pression that the droplet model explains all experimental 
data for different systems. In fact. Figure 4 shows that 
the LDA does not fit the data for the 2DHS density de- 
pendence of dfj,/dp at all. To compare with the data in 
reference [2] we cannot find parameters of a and no which 
would fit the curve even qualitatively. Incidentally, we 
can match the data with an expression which contains a 



It is reasonable to speculate about why the droplet 
model is more appropriate to describe the smaller rg sys- 
tem. It is now a rather common view that the interaction 
can suppress the effect of disorder. As shown in the nu- 
merical simulation, as Ts becomes more important the 
system tends to keep the electron density as uniform as 
possible to minimize the Coulomb interaction. So for the 
strongly interacting heavy mass hole system, it is possible 
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that the droplet state is largely suppressed. 

In summary, we have measured the compressibility of a 
moderately interacting 2D electron system (rg « 0.5-3). 
We revealed qualitative differences between this system 
and that observed earlier in a 2D hole system, where 
interactions are considerably stronger. We found that 
the compressibility data for the 2DES can be described 
in the framework of the droplet state model using the 
local density approximation. 
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grant # DMR 0071969. 
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